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Abstract 

We first review the calculations of low invariant mass dilepton production 
in relativistic heavy ion collisions, using effective hadronic Lagrangians. We 
go through some of the theoretical techniques used in this kinematical region 
and we consider some of the appropriate experimental measurements. Mov¬ 
ing up to the intermediate invariant mass region, we point out some of the 
uncertainties that show up in theoretical estimates. Those originate mainly 
from off-shell effects. Finally, as an application of hadronic chiral Lagrangians 
at finite temperature, we compute the mass shifts and mixing of the oj and cj) 
mesons due to scattering from thermal pions. 


Introduction 

The field of heavy ion collisions is a very active one, straddling high energy and 
nuclear physics. At the upper energy limit of this flourishing area of research, the 
goal is to eventually produce and study a new state of matter in the laboratory: 
the quark-gluon plasma (QGP). That strongly interacting matter in conditions of 
extreme energy densities undergoes a phase transition is in fact a prediction of QCD 
|l|]. To confirm whether the phase transition indeed occurs in relativistic heavy ion 
collisions and that the QGP is formed, one needs a clear signal as a signature for 
the QGP. Many approaches have been suggested to elucidate the existence of this 
elusive state of matter, but unfortunately, no single measurement can be singled- 
out as a “smoking gun” candidate. Instead, it appears that many complementary 
experimental data will require simultaneous analysis |0]. One class of observables 
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that appears especially attractive is that of electromagnetic signals. This owes hrst 
to the fact that such probes essentially suffer little or no hnal state interactions 
and thus constitute reliable carriers able to report on the local conditions at their 
emission site. The calculated emission rates for photons and lepton pairs have been 
shown to strongly depend on the local density and temperature. Those facts were 
established some time ago and several heavy ion experiments specializing in the 
measurement of electromagnetic radiation are either running right now or being 
planned. 

In this talk, we specialize on the topic of dilepton measurements and calculations. 
With respect to photons, the extra degree of freedom associated with a variable 
invariant mass constitutes an advantage that is clear when one considers annihilation 
reactions, for example. With this in mind, a pioneering experiment was that of the 
DLS p], at the disassembled Bevalac at LBNL. Dileptons were also measured by 
the CERES [|| and HELIOS 0 collaborations at the SPS at CERN, and will be 
measured in the PHENIX experiment at RHIC at BNL, and in the HADES 
experiment at GSI 0. 

The emission of lepton pairs from high temperature matter formed in high energy 
nucleus—nucleus collisions is of great theoretical and experimental interest. It can 
signal a change in the properties of hadrons (more precisely, correlation functions) 
in hot hadronic matter as it approaches a chiral symmetry or a quark deconhnement 
phase transition or rapid crossover. Vector mesons are prominent in these studies 
because of their coupling to the electromagnetic current. Most studies have focussed 
on the lighter ones, p, a;, and 0, and on the heavier J/'^. The best signals are those 
that have a characteristic shape or structure. The p meson is very broad in vacuum 
and undoubtedly gets even broader at hnite temperature. The J/V' is narrow; it 
has a whole literature of its own within the held. The uj and 0 mesons are rather 
narrow, but not so narrow that they will all decay after the hot matter has blown 
apart in a high energy collision. This makes them good candidates to study. 


The low mass region 

In general, one can show that the rate for emitting a dilepton pair (or a real photon, 
for that matter) is related to the retarded photon self-energy, at hnite temperature 
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where the photon energy is cu, its three-momentum is k, and T = 1/0. 

A plausible philosophy then consists of performing an order-by-order loop ex¬ 
pansion of the hnite-temperature photon self-energy in the physical helds relevant 
for the problem at hand. One can then establish a correspondence between this 
approach and that consisting of simply using kinetic theory. It can be shown that 
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the difference in the dilepton rate between those two scenarios is indeed small up 
to reasonably high temperatures 0. Systematic evaluation of this sort were carried 
out for the theoretical evaluation of real and virtual photon emission rates from a 
hot gas of mesons [j^, O- We first concentrate on the rates for low mass dileptons 
and recall the assumptions inherent to most “first generation” estimates. We do not 
attempt to do justice to the many calculations in this area: we apologize in advance. 
In the following, we highlight some general aspects. 

Having in mind the energy regime germane to the CERN SPS, a reasonable 
working hypothesis is to first assume that the hadronic medium is meson-dominated. 
This statement is experimentally supported [|T^]. Also, if thermal equilibrium is 
indeed established this would influence the soft sector. This defines a theoretical 
framework where finite-temperature hadronic held theory can be used to evaluate the 
photon self-energy. Its imaginary part can then be related to a dilepton emission 
rate. In practical reality, however, the above-mentioned connection with kinetic 
theory is in fact what is often used. Our “first generation” dilepton emission rate 
then contains the reactions and decays constructed from an ensemble of the lightest 
pseudoscalar (P) and vector mesons {V). Those are (P): tt, t], t]', K, and (R): p, to, (j), 
and K*. The processes included in our rate calculation are then: V{P) 

P-l-P—> 7 *, P-l-R —>- 7 *, and R -|- R —>■ 7 *. The respective sizes of those different 
contribution to the lepton spectrum are compared and discussed in Ref. |TT[ (see 
also Ref. [0). Note that only two-body reactions were first considered, because 
of phase-space arguments. We shall return to this point later. Effective chiral 
Lagrangians were used in this calculation, coupled with Vector Meson Dominance 
[T3[| whenever experimental time-like form factors were unavailable. 

A dynamical model is then needed to translate the rates into actual yields, 
effectively performing a space-time integration. The rates obtained by the procedure 
described above were integrated in a Bjorken model ||T^|, and in an hydrodynamical 


model 


Running the results thus obtained through the detector acceptances. 


one can compare with the measured data. This has been done by several groups 
now, and the results are nicely summarized in Ref. [0. Concentrating on the 


CERES invariant mass coverage, most calculations are seen to underestimate the 
data in the low mass sector (M 300 MeV). One must nevertheless point out 
the remarkable agreement between the different theoretical estimates. The latter 
were obtained within different dynamical frameworks (see and references therein 


for details) and this convergence of theoretical results has been claimed to show 
the insensitivity of the CERES results to dynamical modeling. Going back to our 
rates and recalling that those included many reaction channels (even though a few 
actually dominate), it appears that an independent verification of those estimates 
is desirable. In this context, we focus on two calculations. First, a calculation of 
photon and dilepton production rates from a gas of hot mesons was performed using 
chiral reduction formulae and a virial expansion . The results can be expressed in 


terms of vacuum current-current correlation functions, which can be experimentally 
constrained. Second, using spectral functions extracted from data one can thus 
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Figure 1: A comparison of the dilepton rates at T = 150 MeV calculated using 
experimentally-extracted spectral functions |^] (open circles), chiral reduction formulm 
117] (solid line), and meson reactions and decays (dashed line). 


calculate finite-temperature dilepton emission rates |]^]. A comparison between the 
meson reactions/decays approach outlined earlier and those two calculations appears 
in Fig. 1. It is important to point out that Ref. |]T^ will not have the equivalent of the 
radiative decay channels, as the analysis is done starting from data. Strictly 
speaking, it will then be relevant just below the p pole, where the pion channels 
become large. Also, Ref. [^] does not have the 0 in its formalism: this explains the 
absence of the peak. This being said, those three very different approaches exhibit 
striking agreement. Bear in mind that no arbitrary normalization enters any of 
those calculations. This is not to say that improvements can’t be introduced but 
the unity of the displayed results is nevertheless satisfying. 

Returning to the underestimation of the low mass data by the majority of the 
theoretical efforts, several alternative explanations have so far emerged. Since it is 
not the point of this talk to address this specihc issue with the detail it deserves, we 
simply enumerate some of them. The early suggestion that the CERES data could 
be accounted for by a dropping of the in-medium p mass, has aroused immense 
theoretical interest. This conjecture can be incorporated in a model that hts the data 
. Another line of reasoning consists of trying to make more complete dilepton rate 


assessments. This, combined with a kinetic theory constrained by the available data 
on spectral and multiplicity distributions, appears to be a step towards a resolution 
[0- Note that those rates will have to be reconciled with the 
1. Finally, in-medium baryon effects on the 


of the discrepancy 

alternate approaches described on Fig. 
p spectral function were shown to produce promising results, in connection with the 
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CERES measurements [^. The last two schools of thought do not invoke dropping 
p masses. It is obviously too early to pin down the theoretical explanation which 
will stand the test of time. However, the community is eagerly awaiting the dilepton 
transverse momentum spectra which should help to discriminate between theories. 


The intermediate mass region 


Along with the tantalizing low-mass dilepton excess reported by CERES, an excess 
in the intermediate mass region {m^ < M < has also been seen [^. Before 

attempting to make phenomenological estimates, it is best to exercise caution and 
to tackle this issue purely theoretically. It has been shown in the past that, right 
past its threshold, the reaction vr -|- oi —^ constituted a sizeable fraction of the 

lepton pair spectrum 


22]. This exercise was carried out using an effective chiral 


Lagrangian. 

0 


However, the procedure to derive those is unfortunately not unique 
Now, because we are no longer in the soft sector, our concern is that the 


non-uniqueness of the available chiral Lagrangians will manifest itself, principally 
through off-shell effects. With this in mind, we report here on our findings only 
briefly, leaving a more detailed discussion to be published . 


Here, we therefore carry out a mostly pedagogical exercise: what are the sizes of 
the possible off-shell effects inherent to effective Lagrangians, when used to calculate 
lepton pair production? The definite example we have in mind is again tt + ai —>■ 
e~^e~, evaluated with the help of VMD. We are therefore concerned by the different 
possible manifestations of £,rpai- The source of our concern is illustrated as follows. 
The most general vertex for oi —vrp decay can be written as 


= igi + ig2 -|- igs p''k^ + ig^ + ig^ k^k'’ 


( 2 ) 


where gi{p‘^, k"^) is a form factor, with p^ and k'^ being the the Oi and p four-momenta. 
For on-shell Oi decay, three of the above terms are identically zero, leaving only 
two. The form factors are then chosen such that the experimental width, Ta^^-wpi 
is reproduced. The extra experimental constraint is then the ratio of D to S' wave 
in the final state, which is measured experimentally |^. It is probably fair to say 
that the importance of this hadronic constraint has not quite penetrated the heavy 
ion community, as revealed by a literature survey. When using the above vertex for 
dilepton production, the virtual p is pushed off-shell and one has no direct control 
over such extrapolations. It turns out that for the case at hand, the uncertainty 
actually resides in the form factors p3. We calculate the dilepton rate from the 


reaction above with: (1) an effective chiral Lagrangian constructed by Ko and Rudaz 
based on the SU (2) l x SU {2)r linear a model ; (2) an effective chiral Lagrangian 
where the vector mesons are introduced as massive Yang-Mills fields of the chiral 
symmetry |^; (3) an effective Lagrangian, previously used by Xiong, Shuryak and 
in connection with photon emission rates; (4) an effective Lagrangian 


Brown 27 


developed by Janssen, Holinde, and Speth |^, to address the issue of form factors 
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in the Bonn potential. The dilepton rates evaluated with those interactions are 
displayed in Fig. 2. 

The hrst striking feature of this plot is that the dilepton rates span three or¬ 
ders of magnitude. Off-shell effects are indeed important, as also pointed out by a 

H 


calculation done along similar lines 


For additional help, we can however rely 


on the rates calculated through the experimentally-extracted spectral functions 
used previously. The ti + ai reaction should be a large (if not the largest) con¬ 
tributor to the dilepton signal in the invariant mass region of interest. With the 
available information we can start our examination of the models. Table 1 sum¬ 
marizes their empirical predictions, along with the actually measured quantities. 

0 . 


Parameter set I for the interaction of Ko et al. is taken from Ref. 


Concern¬ 


ing the observables in Table 1, this parameter set does not actually represent the 
lowest parameter set II satishes this criterion. It is interesting to notice that 
a modest variation in for the £t of the on-shell empirical observables translates 
into a large variation of the predicted dilepton rate. Fig. 2 also contains the rate 
deduced from the experimentally-htted spectral functions . Because the process 


discussed here is actually included in the spectral function analysis along with other 
channels (which are known to interfere constructively), the open circles represent a 
strong upper limit for the calculations discussed here. This upper bound is strongly 


I 

> 

CD 

I 

g 


^2 



M (GeV) 


Figure 2: Production rates for intermediate invariant mass dileptons calculated at T = 
150 MeV using the following interactions: (a) and (b) Ko and Rudaz’s chiral Lagrangian 
|]^ (the difference between (a) and (b) is explained in the main text); (c) the effective 
chiral Lagrangian used by Song, Ko, and Gale [22|;(d) Lagrangian of Janssen, Holinde, 
and Speth [^]; (e) Lagrangian used by Xiong, Brown, Shuryak [^]. Also shown is the 
rate calculated from experimentally-fitted spectral functions ||l^ (open circles). 
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Ko et al. 

I 

II 

Song et al. 

Xiong et al. 

Janssen et al. 

DATA 

r(aip7r) 

313.4 

579.1 

400 

400 

400 

~ 400 MeV 

r(ai7r7) 

0.572 

1.171 

2.819 

1.940 

0.312 

0.640 ± 0.246 MeV 

D/S 

0.078 

-0.168 

-0.099 

0.185 

0.045 

-0.09 ±0.03 MeV 


11.2 

9.8 

39.3 

37.3 

11.0 



Table 1: The phenomenology of different models. 


violated by our parameter set II for the Lagrangian of Ref. . To a lesser extent 
this is also the case for the results of calculations with the interaction of Ref. . 
In the remaining models in Table 1, the lowest is obtained with the Lagrangian 
of Ref. Parameter set I of Ko et al. has as good a but the slope of the 


dilepton spectrum departs considerably from that of the spectral function analysis, 
with a violation of the upper bound just below M = 2 GeV. The Lagrangian of 
Xiong et al. produces a high with a much flatter slope. We could deem here 
the interaction of Janssen et al. to be the most satisfactory. It can also be made to 
saturate the spectral function curve with a proper choice of form factors, without 
spoiling its on-shell properties. However, the spectral function analysis can be made 
to be specihc to the ttoi channel and the qualitative aspects discussed here can be 
set on a quantitative footing; we stick here to our pedagogical goal and postpone 
this discussion p4 . 


It is clear which caveats will then await us when we embark on the phenomeno¬ 
logical analysis of the experimental data. The purpose of this exercise was to show 
that using all of the available empirical constraints, one could indeed restore some 
certainty to an apparently ambiguous situation. There are thus reasons to be opti¬ 
mistic about this type of analyses. The experimentally-extracted spectral function 
prove to be helpful, and serve here as baseline calculations. Then, the effective La¬ 
grangian theory can be used to extrapolate into the realm of hnite-temperature and 
many-body effects. 


iu — (j) mixing at finite temperature 

A last application of effective hadronic Lagrangians will cover here the mixing of 
the u and 0 vector mesons. In vacuum the 0 meson is almost entirely ss in its 
valence quark content while the ut meson is almost entirely nonstrange. There is a 
small mixing as evidenced by the observed decay mode 0 —>■ 7rp and by previous 
studies using effective hadronic Lagrangians. We will study the change in masses 
and mixing angles of the a; and 0 at moderate temperatures using only conventional 
ideas. We consider our study a simple extrapolation of known physics which may 
help in deciphering future experiments. Specihcally, we study the scattering of these 


7 























mesons from thermal pions which are the most abundant mesons at temperatures 
below 100 MeV or so. We use effective Lagrangian techniques to model the rela¬ 
tively soft interactions coupling pions, a;, and 0 mesons, htting the parameters to 
known physical quantities. The resulting scattering amplitudes are used in a virial 
expansion to compute the vector meson properties at hnite temperature. 

A survey of the literature on effective hadronic Lagrangians and the Review 
of Particle Physics 0] suggests that the interactions of relevance involve the 3- 
point vertices 0p7r, ojpn, (j)bi7i, and ubin. The vector self—energies are obtained by 
computing one loop diagrams at hnite temperature |^. The interaction involving 


the p meson is described by the Wess-Zumino term in- 


c 






Pp 


■ 77 


'd^ujl + y/2df,ujf, 


( 3 ) 


The octet and singlet helds, ujs and a;^, are expressed in terms of the (vacuum) 
physical helds with a mixing angle Oy- 


Co’s = 0 cos 6v + 00 sin Oy 

OOg = 00 cos 9y — (f) sin 9y (4) 


Ideal mixing is dehned such that the physical 0 meson would not couple to non- 
strange hadrons. This corresponds to 9 y = 6*ideai = tan“^(l/v^) Ri 35.3°. The real 
world is not far from that. Durso [^, for example, hts 39.2°, while the Review 
of Particle Physics quotes 39° based on the Gell-Mann Okubo mass formula. The 
coupling constant g is related to the coupling gyyp used by Gomm, Kaymakcalan, 
and Schechter hy g = —\/ 2 gyyp. Durso hts g'^ = 1.62 ± 0.19 x 10“^ MeV“^ 
with a pseudoscalar mixing angle 9 p = —9.7°. We will insure that our parameters 
reproduce the <p ^ pn decay rate. Accepting Durso’s value of g we ht 9 y = 40.1°. 

The &i(1235) has a branching ratio of more than 50% into oott and less than 
1.5% into 07r. Thus, the bi meson is important for the mass shift of the oo meson at 
hnite temperature, as noticed by Shuryak [^. The interactions are assumed to be 


SU(3) symmetric with SU(3) broken only by mass terms. The Lagrangian is 


r __ uAi I 

^{uj, 4 >)bi 7 T 9 bi * b I 


+ hbj^TT ■ 




Vs 


( 5 ) 


The two coupling constants can be inferred from the decay rate bi — >■ oott and from 
the ratio of the D wave content of the decay amplitude to its S wave content. The 
Review of Particle Physics gives the width 142±8 MeV and the measured D/S ratio 
0.26±0.04. Fitting to central values determines the most likely values of the coupling 
constants to be gt^ = —9.471 GeV and hb^ = 6.642 GeV“^. 

The location of the poles and the mixing angle at hnite temperature are obtained 
by hnding the zeros of the inverse propagator in the (vacuum) physical oo-cj) basis. 

V-\ko, k) = kl-]V-M^ - S(fco, k) 
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is the 2x2 mass matrix at zero temperature. It is diagonal with components 
= m^j and M 22 = The self—energy has contributions from the p meson 
and from the hi meson. The first was calculated by Haglin and Gale |^. The 
second is readily calculated by the usual hnite-temperature rules |^. We restrict 
our attention to the hnite-temperature contribution. We consider only uj and 0 
mesons at rest in the many-body system (k = 0) since that is where temperature 
will have its maximum impact. As |k| increases, many—body effects will decrease, 
and in the limit |k| S> T they will disappear altogether. Finally, the imaginary part 
of the self—energy is small compared to the real part (inclusive of M^) and so we 
do not include it in our present calculation. We hnd 


S(A;o, |k| = 0) = 


cos^ 5v — (cos 5v sin Sy) /2 

— (cos5y sinhy)/2 sin^5y 


(S, + S,J, (7) 


where dy = ^ideai measures the deviation from ideal mixing. The scalar functions 



Figure 3: Mass shifts of the uj and 0 mesons as functions of temperature due to scattering 
from thermal pions. Multiple pion scattering and scattering from other thermal mesons 
should become important above 100 MeV temperature. 


Sp and S 


bi 


are given in Ref. |38 


After diagonalization one hnds the mass shifts as functions of the temperature. 
They are displayed in Fig. 3. Although the results are plotted up to a temperature 
of 150 MeV to see the effect, other scattering processes will come into play above 
100 MeV, as mentioned earlier. The uo mass goes down and a shift of roughly 13 
MeV is reached at a temperature of 140 MeV. The 0 mass increases monotonically 
to about 0.6 MeV above its vacuum value at a temperature of 140 MeV. The mixing 
angles remain uninterestingly small at all temperatures |^|, but the mass shift of 
the UJ is potentially detectable in future experiments. 
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Conclusion 


In conclusion, in our brief description of calculations for the low mass sector we 
have shown agreement between different theoretical approaches to dilepton rate 
calculations. At higher invariant masses, the considerable uncertainties associated 
with off-shell effects have been conhrmed. Nevertheless, we argue that the available 
phenomenology helps considerably in narrowing down the possibilities. Finally, we 
compute the hnite-temperature change in the masses of the lo and 0 mesons and 
their mixing angle. 
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